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' We consider polygonal Wilson loops with null edges in conformal gauge theories. We 

OO \ 

OO . derive an OPE-like expansion when several successive lines of the polygon are becoming 
!>■ ■ 

CN ' aligned. The limit corresponds to a coUinear, or multicoUinear, limit and we explain the 
^ ! systematics of all the subleading corrections, going beyond the leading terms that were 
5 ■ previously considered. These subleading corrections are governed by excitations of high 
spin operators, or excitations of a flux tube that goes between two Wilson lines. The 
^ . discussion is valid for any conformal gauge theory, for any coupling and in any dimension. 
■ For M = A super Yang Mills we check this expansion at strong coupling and at two 

loops at weak coupling . We also make predictions for the remainder function at higher 
loops. 

In the process, we also derived a new version for the TBA integral equations that 
determine the strong coupling answer and present the area as the associated Yang- Yang 
functional. 
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1. Introduction 

The Operator Product Expansion (OPE) is a powerful tool for studying correlation 
functions in conformal field theories. The expansion is controlled by the spectrum of local 
operators of the theory. In particular, their dimensions control the powers of the expansion 
parameter in the OPE, which is simply the separation between operators. In some cases, 
such as two dimensional minimal models, it is possible to completely fix the correlation 
functions by demanding this property in all possible channels [0. 

In this paper we derive a similar OPE-like expansion for polygonal Wilson loops with 
light-like edges. Our main motivation is the relation between Wilson loops and amplitudes 
in = 4 super Yang Mills . Furthermore, similar looking Wilson loops appear 
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in a variety of high energy processes in gauge theories. The OPE expansion is vahd for any 
conformal field theory and in any dimension where we can define nuU polygonal Wilson 
loops. It is valid whenever the Wilson lines produce a conserved flux that cannot be 
screened. It is useful to consider flrst large gauge theories in the planar approximation. 
We will later make some remarks beyond the planar limit. One challenging aspect of the 
problem is that Wilson loops with light-like edges are eminently Lorentzian observables 
without an obvious Euclidean counterpart. 



Fig. 1: We consider a general Wilson loop with null edges. We select two of its 
edges and we place and imaginary cut along a line connecting these two selected 
edges. We then expand the answer in terms of states propagating across this cut. 

Our OPE expansion is performed as follows, see fig. 1. First we select two non- 
consecutive lines of the polygonal Wilson loop. We then cut the Wilson loop into a top 
part and a bottom part and we expand it in terms of the states propagating through the 
cut. The states that propagate through the cut correspond to excitations of a flux tube 
that ends on two light-like lines. Fortunately these states have been considered before. In 
fact, they are excitations around the infinite spin limit of high spin operators. In theories 
with gravity /string duals they are excitations around the spinning string (GKP) considered 
in . In A/" = 4 super Yang Mills one can compute these dimensions for all values of the 
coupling [^,|^. Thus, we find that the Wilson loop has an OPE-like expansion in terms 
of operators, or states, which are excitations of the infinite spin limit of the GKP string. 
These states can also be viewed as created by local operator insertions along a null Wilson 
line. 
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Fig. 2: The collinear expansion including subleading terms. The first term is the 
usual statement that in the collinear limit we recover the Wilson loop with one 
less line. The second term corresponds to the insertion of one operator along the 
contour, the second to two operators, etc. 



The simplest example of this OPE is a collinear limit, where two consecutive lines 
become parallel, see fig. 2. This collinear limit was considered previously and its leading 
divergent and constant terms were understood ITT]. These leading terms are determined 
by the special conformal symmetry anomalous Ward identities We are now saying 

that we can continue the expansion into the subleading terms. The existence of this 
expansion places constraints on the possible form of the Wilson loop correlator. In fact, 
it constrains the so called remainder function [^|TB,n,3| which contains the conformal 



invariant information on the Wilson loop correlator. We will demonstrate this explicitly 
below for the case of the six sided Wilson loop in A/" = 4 super Yang Mills. In fact we will 
verify that the strong coupling answer and the two loop weak coupling answer have the 
form required by the OPE. We will also make a prediction, using for the expansion 
for all values of the coupling. At strong coupling we can go further and check that the 
expansion has the required form for general polygons. 

It is quite likely, that the existence of this expansion in all possible channels, to- 
gether with integrability would completely constrain the correlator. This would simply 
be the extension of the Bootstrap procedure for correlation functions to Wilson loop 
correlators. 

Our paper is organized as follows. In section 2 we discuss the symmetries of the 
problem. We will use them heavily for deriving the existence of the expansion. The key 
idea is to identify a "Hamiltonian" which organizes the expansion. In section 3 we explain 
the form of the OPE expansion for Wilson loops. In section 4 we discuss the form of the 
expansion for the case of the hexagonal Wilson loop at strong and weak coupling. Finally 
in section 5 we briefiy discuss higher order predictions. 



^ Though this statement is believed to be true, to our knowledge it hasn't been proven. We 
will see that it follows from our discussion at the end of section 3. 
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The paper contains several appendices. In appendix A we describe our setup in terms 
of momentum twistors. In appendix B we present further arguments for the identification 
of the states we encounter in the expansion and excitations around the flux tube connecting 
two null Wilson lines. In appendix C we discuss the dispersion relation of the eigenvalues 
of the "Hamiltonian" . In appendix D we present some details on the one loop computation 
of some "form factors" . 

Finally, in appendices E and F we give a new form for the integral equations described 
in [ITSlJT^] . In the new form, the equations involve only physical cross ratios. Furthermore, 
the expression for the area is given by the Yang- Yang functional of the associated TBA-like 
integral equations. This form is particularly useful for analyzing our limit. It is very likely 
that it will also be useful for further studies on integrability and Wilson loops/amplitudes. 
Thus these appendices can be read on their own. 

2. Symmetries 

The expansion we are studying is governed by the symmetries preserved by two null 
lines and also by the symmetries preserved by the square Wilson loop, i.e. a Wilson loop 
with four null sides. 

To justify the need to understand the symmetries, let us recall first how symmetries 
constrain and determine the ordinary OPE for correlation functions. In that case the main 
symmetry in question is dilatations. If we have a correlator {0{1)0{2) ...0{n)) . Then we 
can set 0{1) at the origin and consider the limit where 0{2) approaches it. This operation 
corresponds to a dilatation applied to operators 0(1) and 0(2), but not to the rest of 
the operators. This becomes more clear if we map this to the (Euclidean) cylinder. Then 
0(1) is at Ti = — oo and 0(2) is at some value of T2 and the rest are at fixed values 
of Ti. Here r is the (Euclidean) time coordinate on the cylinder. We also demand that 
T2 < mm{Ti}. The OPE expansion follows by cutting the cylinder between T2 and min{ri} 
and inserting a complete basis of eigenstates of the dilatation operator. The fact that the 
OPE is convergent follows from the asymptotic growth of statesil. 

^ For a d dimensional CFT we have logA'^(A) (const ant) A <c A as A — > 00. 
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2.1. The symmetries preserved by null lines and the square Wilson loop 

Let us start with two generic non-intersecting null lines and consider the symmetries 
that leave these two null lines invariant. Working in R^'^, by conformal transformations we 
send a null line to null infinity and the other to a null line extended along x~ and passing 
through the origin. We can do this so that both null lines are in axi R^'^ subspace, with one 
of the lines at null infinity. The symmetries that leave this configuration invariant include 
dilatations D, boosts in x^, M"' , and a rotation in the two transverse directions M^^. 
In addition we have translations P_ along x~ and special conformal transformations . 
The symmetries P_, K~ and D + M+~ form an SL{2,R) subgroup, which commutes 
with the two remaining abelian symmetries: D — M+~ and M^^. Note that these two 
remaining symmetries leave individual points on the lines invariant, while the SL{2, R) 
symmetry maps points on a given line to other points on the same line, but does not leave 
them invariant in a pointwise fashion. 

For this counting of symmetries it is important that we have the line at null infinity. If 
we only had the line that goes through the origin we would have more symmetry generators 
such as K^, K'^, and boosts M+* ~ x'^di — x^dr^- , where z is a transverse index. However, 
these transformations do not leave the line at infinity invariant. This can be checked by 
using an inversion to map the line at infinity to the origin. For example, this turns into 
which indeed does not leave the line at the origin fixed. These symmetries preserve one 
of the null lines but move the other one. So they are not true symmetries of the two line 
system. Nevertheless we will see that they have some interesting consequences. 



Fig. 3: We start with two null lines, denoted by the solid lines. By picking a 
point A on the first line, we determine a unique point C on the second line which 
is light-hke separated from A. Similarly for B. In this way we construct a square. 
In the figure we projected onto two spatial directions and suppressed the time 
direction. 
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Fig. 4: (a) The square Wilson loop projected onto two spatial directions. We 
suppressed the time direction. The "Hamiltonian" symmetry dr moves points 
along the arrow direction. It leaves points on the thin black lines fixed, but moves 
points along the thick red lines. The "momentum" symmetry da is also indicated, 
(b) We inserted an operator on the bottom line. We can integrate the operator 
along the bottom edge in such a way that it has a definite momentum with respect 
to the da symmetry. If it creates a single particle state, then it will also have a 
definite energy, (c) We have mapped the square to an i?^'^ subspace by a conformal 
transformation that sent a cusp to spatial infinity. Two of the null lines are at null 
infinity. One of the cusps is at the origin. Operators inserted along the top and 
bottom lines in (b) correspond to operator insertions which are spacelike separated, 
and indicated by a cross. The "Hamiltonian" and "momentum" correspond to the 
generators D it . 

Once we have understood the symmetries of two lines, we now want to understand 
the symmetries of the square. To define the square we need to pick two points along one 
of the null lines where the vertices of the square will be sitting. Once we pick a point 
along a null line to be a vertex of the square, the corresponding point along the other null 
line is automatically fixed, since it is the point on the other null line that is null separated 
with respect to the point chosen on the first line, see fig. 3. Thus, picking a square, now 
corresponds to picking two points along a null line. This breaks the SL(2, R) symmetry 
group to a single generator. We are going to view this remaining generator of SL{2, R) 
as a "Hamiltonian" . This Hamiltonian moves points on both segments of the original null 
lines towards the same side of the square. In fact, we will see that, because the new sides 
of the square are spacelike separated, this "Hamiltonian" is a Euclidean Hamiltonian such 
that the evolution operator looks like e""^^, as opposed to e"*"^^, where E are the real 
energy eigenvalues. The second non-compact symmetry moves points along the bottom 
and top lines of the square, see fig. 4(a). We can view this as a "momentum" generator. 
In fact, the "Hamiltonian" has a continuous spectrum due to the existence of this second 
non-compact symmetry. 
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In summary, the square has three commuting symmetries. Two noncompact symme- 
tries and one transverse rotation symmetry. In (2,2) signature the transverse symmetry 
becomes a boost and is also non-compact. 

The symmetries of the square can also be understood directly as follows. By a confor- 
mal transformation we can map three of the points to the boundary of Minkwoski space. 
One of the points is mapped to spatial infinity, and the other two are mapped to some 
points on null infinity. The fourth point can be set at the origin. In that case we just 
have a null cusp with lines along x^x~ = 0, > 0. See fig. 4(c). Let us understand 
the symmetries that leave this invariant. We clearly have three commuting symmetries: 
dilatations, boosts in the plane and the rotation in the transverse plane. In fact, these 
are all the symmetries of this Wilson loop. It is fairly clear that there are no further 
Poincare symmetries of R^'^ that leave it invariant. There are no special conformal sym- 
metries that leave it invariant because any special conformal symmetry moves the point 
at spatial infinity, which is one of the vertices of the polygon. 
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Table 1: Symmetries preserved by one and two lines and the square. 
In the last line we introduced some notation for each of the three symmetries. 

It is also useful to understand how the "Hamiltonian" acts on the "bottom" side of 
the square, see fig. 4. It leaves it invariant on a point by point basis and it rescales the 
transverse directions. More precisely, if we say that this bottom side is along x^ , then 
it rescales x~ — )■ X^x~ and y± — )■ Xy±, see fig. 4(c). This action is sometimes called 
"twist" operator. More concretely, we will find that in our expansion procedure we will 
need to consider operators inserted along this "bottom" null line, see fig. 4(b) (c). Similar 



operators arise also in the study of high energy processes in QCD, for a review see ||17|| . 
Such an operator is characterized by a "momentum" along the direction of the line. This 
is the same "momentum" generator that we mentioned above. In the QCD literature on 
operators involving null Wilson lines this generator is called the spin of the state (e.g. in 



||17|| ). It turns out that if we insert a local operator along the null Wilson line, this operator 
will not be an eigenstate of the twist operator (or our "Hamiltonian" ) . However, single 
particle states with definite "momentum" are indeed eigenstates. 

The states that are contributing to the OPE are the states created by inserting local 
operators along the bottom side of the Wilson loop and propagating all the way to the top, 
see fig. 4(b). These states have a continuous spectrum, due to the existence of a second 
non-compact symmetry. This is in contrast to the ordinary operator product expansion 
where the spectrum of dimensions is discrete. However, due to the momentum symmetry 
we can certainly isolate the contributions from "single particle" states propagating along 
the square. 

3. The collinear operator product expansion 

In this section we state more precisely how to perform the OPE expansion. In the case 
of ordinary correlation functions the expansion can be organized in terms of the dimensions 
and the spins of the intermediate operators. As we mentioned above, the operators that 
appear in our expansion are characterized by their dimension (or energy) and also by 
another continuous label which we will call the momentum. There is a third discrete label 
which is the spin in the transverse dimension. In order to isolate the various quantum 
numbers it is convenient to introduce the following three parameter family of polygons. 

3.1. A three parameter family of polygons 




(a) (b) (c) 



Fig. 5: (a) Given a general polygon we select two segments. We prolong them 
and choose two points on one of these segments A and D. These determine points 
B and C on the other line. We use them to form a reference square ABCD. We 
now cut the polygon into a top part and a bottom part. We apply a combination 
of symmetries of the square to the bottom part. In the limit r — >■ oo the bottom 
part is flattened out into the straight edge of the reference square. In that limit 
we get the polygon in (b) which we will call the "top" polygon. Alternatively, we 
could have applied the symmetry to the top half of the polygon in (a). In the limit 
we would get the polygon in (c) which we call the "bottom" polygon. 
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We would now like to introduce an interesting family of polygons. The family is spec- 
ified as follows. Let us assume that all the points on the polygon are spacelike separated, 
except for the obvious null separations that define the polygon. We first pick two non- 
consecutive null sides of the polygon, and extend them to two null lines. We then pick two 
points along the null line which lie outside the segment that belongs to the polygon on 
that null line. These two points, together with the two null lines define a certain reference 
square. The symmetries of the square involve the three symmetries we discussed above, 
see fig. 5. Let us say the two lines we picked are lines i and j. This splits the polygon 
into a "top" part and a "bottom" part. We now consider a symmetry generator M which 
leaves the reference square invariant. We leave the "top" invariant and act with M on the 
bottom part of the original polygon. The cross-ratios for this family of polygons can be 
described more explicitly in terms of momentum twistors, see appendix A. 

The family has three parameters M = e^^-+'^^-+'^^^ = Q-rE+iaP+i^s _ r^^ie three 
parameters r, a, are conjugate to the three symmetries of the problem. The r parameter 
is the one that is coupling to the symmetry that was included in the SL{2, R) symmetry of 
the two null lines. When we send r — > oo we are mapping the whole set of points between 
lines j and i to one of the sides of the square. We can view this as a multicoUinear limit, see 
fig. 5. The coUinear limit corresponds to the special case where we have only two segments 
between segments j and i and we take the r — > oo limit, see fig. 2. 

The idea is that we can expand the whole Wilson loop in terms of eigenstates of these 
three operators. These eigenstates corresponds to excitations of the flux tube (or string) 
that goes between the two selected null lines. Thus, we are cutting the Wilson loop. The 
bottom and top parts of the loop correspond to two particular superpositions of states. 
Introducing the family of polygons is helpful for isolating the contributions from different 
eigenstates. To perform this expansion we are simply using the conformal symmetry of 
the problem. Thus, the expansion makes sense for null polygons in any dimension where 
we can define Wilson loops and for any value of the coupling. 

Again, given a polygon we fix a reference square and introduce a family of polygons 
parametrized by r, a and 0. These coordinates are associated to the unbroken symmetries 
of the square. The original polygon may be associated to the point t = a = (f) = but that 
is of course not important, we then expand at large r and resum to get r = if we want. 
The choice of reference square is the analog of the choice of point for an operator product 
expansion and the points r, a and are the analog of the coordinates with respect to that 
chosen point. 
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So far we are describing the kinematics of the expansion we want to define. In order 
to flesh this out a bit more we wiU need to understand some more aspects of the dynamics 
so that we have a clearer picture of the states that appear in the expansion. 



3.2. States propagating in the expansion 

What is interesting about this family of polygons is that we can view the symmetry 
conjugate to r as a Hamiltonian. In that case we expect that the result has an expression 
of the form 



where n is some set of labels for the different "energy" eigenstates. We have also em- 
phasized the fact that the coefficients Cn factorize into two contributions which are the 
overlap of the intermediate state with the top or bottom polygon. In order to make a 
precise statement we will have to take into account that UV divergencies of the Wilson 
loop break these symmetries. Fortunately, this breaking is well understood and we will be 
able to take care of the those effects in a simple fashion. So, let us treat it as an exact 
symmetry for the time being. 

Now, let us explain better what kind of states we expect in (|3.1|) . In the case of the 
ordinary OPE of local operators we can surround the two operators we are considering by 
a three sphere. We can then view the states appearing in the OPE as energy eigenstates 
of the theory on the sphere. To be more precise, we could consider the conformal field 
theory in R x and find the energy eigenstates. 

In our case, the procedure is very similar, we could consider a constant r surface. 
We see that the surface is pierced by two null Wilson lines, see fig. 1. The Wilson lines 
that pierce the surface create a color electric fiux tube. The states can be understood 
as excitations of this fiux tube. In other words, the fiux with no excitations gives a fiux 
vacuum. We need to understand the excitations of the field theory around this vacuum. 

More explicitly the states propagating can be understood as follows. We start with 
the Yang Mills theory on i? x 5""^. We now add two null Wilson lines moving along a great 
circle of the S^, see fig. 6(a). If (f is the coordinate of that circle, then the lines are along 
t = (f and t = (f + n where t is the time coordinate. The generator A — S = idt + id^p 
is a symmetry of the configuration and we will call it the "Energy" . It is also called the 
"twist" . One line corresponds to a quark and the other to an anti-quark. These lines create 
a color electric field. The configuration is invariant under the "momentum" symmetry, see 





im 



(3.1) 
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(a) (b) 

Fig. 6: (a) Field theory on R x in the presence of two null Wilson lines. 
Dotted lines are in the back of the cylinder. Of course the circle of the cylinder is 
really an S^. (b) Double analytic continuation sends the initial and finial states of 
(a) to two null lines, the two blue null lines. Different states will contain different 
insertions along these two blue null lines. 



||18|| , generated by P = —ida- This is not an obvious geometric symmetry of R x S^, it 
is a combination involving conformal Killing vectors. The flux is extended along the a 
direction and has a constant energy density along a. Thus its energy diverges. The energy 



density is simply the cusp anomalous dimension [|T^. The flux emanates from a point in 
the extra S'^ C and the transverse SO (2) symmetry is the rotation that leaves this point 
invariant. In appendix B we show how to choose a conformal frame where both the r and 
a translations are manifest. This is the ground state of the conflgurations we consider. 
We can then add excitations which propagate on top of this conflguration. In the planar 
theory these are excitations of the color electric flux and the indices of the particles we 
create are contracted with the indices of the background flux. This is a Lorentzian picture 
for the states appearing in the OPE. See appendix B for further discussion. 

We can perform an analytic continuation of the conflguration with two Wilson lines 
on i? X S"^ to a conflguration with four Wilson lines describing the square Wilson loop, 
see flg. 6(b). Upon this transformation the Energy generator becomes a generator which 
acts in a Euclidean fashion. It becomes the dr generator we mentioned before for the 
square. The two lines of the Lorentzian picture become two opposite lines of the square. 
The excitations are produced by adding operators on the other two lines of the square, see 
flg. 6(b). 
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Fig. 7: Expansion of a Wilson loop in terms of states propagating on the square. 
The first term corresponds to the expectation values of the top and bottom Wilson 
loops. We can define the expectation value of the square as being one. The second 
term corresponds to the exchange of a single particle. The top and bottom Wilson 
loops give rise to the factors C*°p and C''°"°" in the OPE expansion. The third 
term contains two particle states, and possible bound states, etc. 



These states form representations under the two other commuting symmetries. One 
of them, da, is a non-compact symmetry. Thus, we expect a continuous spectrum of 
excitations. We can thus separate states according to their "momentum" in the a direction. 
In addition, we have some angular momentum number m under the transverse rotation. 
Thus a more refined statement has the form 

{W) = j dne-^^"+'P"''+''^"'^Cr, (3.2) 

Now that we have defined this family we can simply take the large r limit. This would 
select the lowest lying states from (|3.2|) . But the expansion makes sense for any r and we 
conjecture that it converges for all values of r such that all points in the polygon remain 
spacelike separated for all a and (p for that fixed r. 

The first subleading correction comes from a single excitation that is moving on the 
background of the square Wilson loop. This corresponds to the lowest lying excitation 
of the color electric fiux we discussed above. In other words, taking the leading order 
limit corresponds to replacing the whole bottom of the Wilson loop by the single line 
of the square, fig. 5(b). Taking into account the fact that the line is not really straight 
corresponds to insertions of F^i along that line, where + is the direction along the line. This 
is the case, because deformations of a Wilson loop can be understood in terms of insertions 
of such operators along the original loop. Thus, the simplest deformation corresponds to 
a single insertion of the field strength F^i along the loop. By an insertion along the loop, 
what we really mean is that we have an expression of the form 

Tr[Pef^^F+,{x)] (3.3) 
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where F is inside the trace and inserted along the contour at point x. ' denotes the 

integral of the connection along the loop starting at x, going around the loop and ending 
again at x. 

Note that F is inserted along a null line. We can now consider the generator A — S — 
D — M^ which leaves points on this line invariant. This is the so called "twist" generator. 
In fact, operators which correspond to Wilson lines along a null direction with operators 
inserted on them were analyzed quite extensively in gauge theories, including QCD, because 
they govern many interesting high energy processes, going back to the classic analysis of 
deep inelastic scattering. See IJl^ for a review. 



Thus, one of the interesting points is that the operators that appear are fairly well 
understood and have been studied in the past. Note that the operator (|3.3|) breaks the 
cr translation symmetry, since moves points along the null line where Fj^i is inserted. 
Thus, we can consider superpositions with definite momentum p along the a direction. In 
addition, it has charge ±1 under the transverse SO{2) symmetry. There is a unique state 
with momentum p and thus this state is automatically an eigenstate of the Hamiltonian 
dr, with energy t{p). We can think of e{p) as the single particle dispersion relation. Our 
discussion so far, has been valid for any conformal planar gauge theory. In the particular 
case of A/" = 4 super Yang mills the function e(p. A) has been computed exactly in 0. The 
square has a symmetry under the exchange of r and cr. This amounts to a "Wick" rotation 
which exchanges the spatial and time direction. This implies certain constraints on the 
dispersion relation. More can be found in appendix C. 

For general non-planar conformal theories this discussion continues to be valid as long 
as the flux is conserved. For example, in non planar M = A super Yang Mills the discussion 
continues to be valid if we consider Wilson lines in the fundamental. The flux vacuum is 
well deflned, and so are its excitations. If we consider Wilson lines in the adjoint, then the 
flux can be screened and the discussion will need some modiflcations. 

In summary, the leading order expression in the large r limit is given byi 

\og{W) = log(Vr)*°P + log(W^)^°"°'" + 5Z / ^P^^ A)e±^'^e*P"e-^(P'^)" + ■ ■ ■ (3.4) 

where C±(p, A) are some unknown functions, analogous to OPE coefficients. And 
log(VF)*°P' are the expectation values of the top and bottom Wilson loops, see fig. 5 



^ Taking the logarithm reorganizes the expansion in the usual fashion. 
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(b)(c). The function e{p, A) is the exact dimension (or twist) for a single insertion of F^i 
along a null line. In fact, C± = C^"^ x (7bottom ^j^g];g q^p g^^id (7bottom correspond to the 
expectation values of the Wilson loop contours in fig. 7 with an F^i insertion on the top or 
bottom ImM. In other words, as usual, C^"^ is the overlap between the intermediate state 
that is propagating and the state created by top part of the Wilson loop while (7bottom jg 
the overlap with the bottom Wilson loop. 

3. 3. Taking care of violations of conformal symmetry: the remainder function 

So far we have ignored the UV divergencies of the Wilson loop and we have treated 
the symmetries as if they were unbroken. Since the breaking of the symmetries is a well 
understood phenomenon |Tl[], it is clear that we should be able to take this breaking 



into account. In other words, there is an anomalous ward identity which tell us how the 
Wilson loop changes when we apply a special conformal transformation. Thus, we could go 
through the above argument adding the corresponding terms due to the anomalous ward 
identity. 

We will follow a slightly different route which we found more convenient. This is based 
on the observation that the anomalous ward identity is also obeyed by the same Wilson 
loop correlator but in a free U{1) theory. More precisely, we simply need to take the free 
U{1) result and replace the cusp anomalous dimension of the U(l) theory, Fi ^cusp, by the 
full cusp anomalous dimension, Vcusp of the interacting theory in question. So we can 
write 

^ cusp 

(^)^(i) = [(^)u{i)\ = e^^-^^d) (3.5) 

where wu{i) is the result for a U{1) theory. This function is also related to the one loop 
maximally helicity violating amplitude once the tree level contribution is stripped out 
|l^,|10[- It is given by a single gluon (or rather photon) exchange between all pairs of lines. 
It is a completely explicit function of the distances between various cusps and we will not 
need its explicit expressioni. The tilde in ^7(1) just reminds us that we have put the cusp 
anomalous dimension of the interacting theory. 



^ More precisely, they are the ratio of the expectation value with an insertion divided by the 
expectation value without the insertion. 



^ Its explicit form can be found in formula (4.58) of [10|. In this context this function is known 
as "the BDS expression". 
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Thus the ratio 

{W) 



R = loe 



(3.6) 



is a conformal invariant function. This is also caUed the remainder function. By definition, 
its first non-zero contribution is at two loops [p!3| ,|6|,p^ . This definition clearly takes care 
of the divergence problem and leads to an explicitly conformal invariant answer. Though, 
this ratio takes care of the double logarithmic divergencies, there are single logarithmic 
divergencies that should be taken care of. These are regulator dependent. Thus, we adjust 
the regulator in the U{1) theory so that the single logarithmic divergencies match those of 
the full theory. Similarly there is a finite constant which is proportional to the number of 
cusps and is also regulator dependent i. These single log divergent terms or the constant 
terms are not important for this paper. 

This ratio gives a nice conformal invariant expression, however, we should understand 
how it modifies the expectations from the point of view of the Hamiltonian interpretation 
of the family of polygons and the expansion for large r. 

For this purpose, it is useful to note that does have its own Hamiltonian 

interpretation since it is a computation in a U{1) theory. Thus, we have an expansion of 
the form in ( |3.2| ) where the anomalous dimensions vanish, and the energies appearing in 
( p.2| ) are simply the twist of the corresponding operators in the free U{1) theory. In the 
free U{1) theory the twists are just integers. Thus the expansion in (|3.6| ) contains the ratio 
of the two expansions, or difference once we take the log. 

In other words, the final version of the expansion takes the form 



(3.7) 



where the subindex indicates that we are considering the Uil) theory. and i^^ottom 
are the remainder functions for the two polygons in figure fig. 5(b),(c). We are ignoring 
here a possible constant piece which is independent of the kinematics. 

Now, in practice, we can distinguish the terms that come from the free theory because 
the energies are those of a free theory. 



^ We thank E. Sokatchev for pointing out the omission of these single logarithmic terms from 
the first version of this paper. 
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For example, if we consider the terms arising from a single insertion of the field F^i 
then the expansion takes the simple form 



R = i^top ^ ^bottom ^ J2 e^'^ j 



dpC±{p, A)e" 



-e(p,A)T+ipcr _ g — T 



dpClip)e'P'' 



(3. 



The last term is the contribution from the U{1) theory and it has a simple r dependence. 
The first term includes the full dispersion relation and this leads to a more complicated r 
dependence. In the next section we will check that for the case of an hexagon the strong 

Tl| has the form 



coupling answer |]T5| , as well as the two loop weak coupling answer []T^ 
predicted by ( p.8|) . We will also discuss the general structure of n-sided polygons at strong 
coupling. We see that (|3.8| ) implies a particular structure for the Wilson loop correlator. 

The dispersion relation e(p, A) in (|3.8| ) has been computed for all values of the coupling 
in 0. Using that result, then ( |3.8[ ) gives a prediction for the Wilson loop (or the MHV 
amplitude if they are equal) that should hold for all values of the coupling. Notice that 
we do not know what C is for all values of A. Nevertheless we have a concrete prediction 
in terms of the r and a dependence for the answer. 

In appendix C we make further remarks on the spectrum of excitations around the 
flux tube. It would be useful for the reader to consult this appendix if he/she is not familiar 
with these excitations. 



3.4- Taking care of violations of conformal symmetry without using the U{1) or one loop 



answer 



It would be nice to be able to strip off completely the ^7(1) theory contribution from 
our expansion. This is actually possible, by regulating the anomalies in a different way. 
It is useful first to understand how this can be done in a U{1) theory. We start with the 
following ratio of Wilson loop expectation values in a U{1) theory. 



ruii) = log 



{W)u(l) {Wsquare)u{l) 



(3.9) 



is finite and conformal invariant, if 1^*°p is the Wilson loop with the bottom part replaced 
by the single bottom line of the square and VF*^"**"™ is the Wilson loop with the top part 
replaced by the single top line of the square, see fig. 5(b),(c). 
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(. <Wsquare> 
<Wbottoin> 



<Wsquare> 
<Wtop> 



(a) (b) 




(c) (d) 



Fig. 8: (a) The original contour with the reference square given by its vertices 
ABCD. (b) The OPE expansion in the U{1) theory can be written in terms of 
r[/(i) which is computed by a single photon exchange between the top and bottom 
Wilson loops in this figure. In (c) and (d) we see the same, but with a choice 
of points that removes the divergencies. Points D and B of the reference square 
coincide with points on the original polygon. 



The ^7(1) Wilson loop is the double integral of the free propagator 

\og{W)uii) = \ i dx^ I dx"'G^^{x,x') (3.10) 
Jw Jw 

where we regularized the integral by shifting W slightly in the second contour integral. 
Then it is easy to see that 



ru{i)=f dx^ i dx'^Gf.^ix.x') (3.11) 

Here the contours Wsquare — Wtop and Wgquare — Wbottom are well separated and one is 
tempted to say that this expression is finite. However, there is a remaining divergence from 
propagators that connect the short portion of the null line in the bottom diagram with the 
segments that touch the extension of same null line in the top , see fig. 8. This leads to 
divergent terms that depend only on the location of the points along the null lines though 
AD , for examplell. Such divergencies are independent of a. They can be interpreted as 



^ There is a single log divergence proportional to the log of the cross ratio of the four points that 
are sitting along the null line passing through AD. In other words, we have a term log fi log , 

AD AD 

where AD are points on the reference square and A, D are points on the original polygon along 
the same line, see fig. 8(a). There a similar one from BC. These cross ratios depend only on r 
and not on a. They vanish if D = I? as in fig. 8(d). 
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coming from a condensate of Goldstone bosons of the broken SL(2, R) symmetry that acts 
on r. These divergencies can also be removed by choosing points B and D (or A and C) of 
the reference polygon to coincide with the vertices of the original polygon, see fig. 8(c) (d). 
Thus, we can either remove these divergencies by an appropriate choice of the reference 
polygon, or we can just simply ignore them because they are independent of a. 
The expansion of rj/^i) has the correct form to be identified with 

ruii) = j dnC^e-^'^+'P>+'<'^ (3.12) 

Indeed we can expand the propagator in eigenfunctions of the symmetries of the square, 
including the action of M on the bottom Wilson loop: 



G ^i,{x,x'^dx^ dx'^ = I dn dx^-^^{x,n) dx'^'ip^i^x' ,n) e 
and compute the sources 



c^°P{n)=i dx^'^^ix.n) (3.13) 

^ ^square ~^top 

and c'l^ttomi''^) so ^^^^ '^n = ctop{n)cbottom{n). In the next section we check this result 
for the hexagon Wilson loop. Note that this gives an explicit formula for the coefficients 
qO _ ^cuap ^0 appearing in the expansion (13.81). 

-L l,cusp ' ^ 

With this ^7(1) result as an inspiration we are lead to define an alternative conformal 
invariant and finite expression! 

{W) {Wsquare) 



log 



(3.14) 



ottom I J 

where now all expectation values are in the full interacting theory. This leads to a simple 
expansion involving only the physical excitations 

dn c^e-^"^+^P"''+^"""'^ . (3.15) 



The combination in ( p.l4|) is subtracting the contribution from the unexcited flux tube, 
or flux vacuum contribution. The remainder function deflned in ( |3.6| ) is simply given by 
R = Ri R2 + r — ^ij(^iy Notice that the fact that ( p.l4|) has an expansion in terms of 



subleading terms only also explains why the leading terms in the ordinary coUinear limit 
are flxed by conformal symmetry. In that case VFbottom is a pentagon whose expectation 
value is flxed by the anomalous ward identities [0. This pentagon leads to the terms in 
the splitting function computed previously [1^ . 



^ Again we choose D = D and B = B to eliminate a single log divergence. Alternatively, we 
could consider a more general reference square and ignore the a independent divergence. 
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4. Detailed checks for the hexagon 

A simple case that we can consider is the hexagonal Wilson loop. In this case there 
are three cross ratios ui,U2,us- Since our family involves three parameters, then it is clear 
that we can parametrize the three cross ratios in terms of the three parameters r, a, (p. In 
this case, the limit r — t- oo is a coUinear limit which leaves behind a pentagon, which does 
not have any remaining cross ratios. 

There are different ways to parametrize the Ui. An important choice is the choice 
of two opposite lines of the hexagon. There are three different ways of doing this which 
represent the three different channels in term of which we can expand the correlator. Once 
we chose a pair of opposite lines we can now choose the reference square in slightly different 
ways which amount to the action of the generators of SL{2, R) that leaves the two lines 
fixed. Different choices give slightly different parametrizations which only affect results at 
higher orders. These different choices change the answer in a predictable way, as is the 
case for the ordinary OPE. We can make a choice so that the cross ratios are, see appendix 
F.3, 

1 

U2 = 7^ 



cosh^ T 

e'^ sinh r tanh r 
^ 2(— cos ^ + cosh r cosh (j) ('^■-'-) 
e~'^ sinh r tanh r 

^3 



2(— cos (p + cosh r cosh a) 

The expression for ui and us appears rather complicated, but one may note that ui/us = 
e^" and that the parameter \x appearing in the integral equations in [jl5| is /i = —e^'^. In 
(2,2) signature we can set [i = and cos^ — )■ — cosh/. 

Once we have this parametrization we can expand the remainder function for large r. 

J^.l. Expansion of the hexagon at strong coupling 

In this section we quote the result of the expansion of the hexagon at strong coupling. 
We leave the details to appendix F, where we show that the remainder function at strong 
coupling does indeed admit an expansion of the form (p.8D involving a reasonable spectrum 



of operators, thanks to the cancellation of several unwanted terms. 
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We expand the answer to order e We find 



R 



Ri + R^ + R2 + • • • 

— cos (/)e~'^ (cosher log [2 cosh a] — asinha) 



Ri 



^V2 



R2 




2tt [sinh(20 + zO)]2 




(4.2) 



Let us explain the interpretation of the various terms. The terms multiplied by cos(/) 
correspond to the propagation of a particle created by the insertion of an excitation F+i 
which carries unit charge under the transverse SO (2). The Ri term has energy (or twist) 
one, independent of the momentum. This is then interpreted as coming from the U{1) 



sponding particle at strong coupling. At strong coupling this particle has a relativistic 
dispersion relation with mass \/2; its momentum is p = y^sinh^. The term R2 contains 
terms going like e"^"^ which come from the contribution of the U{1) theory. One of them 
is independent of (p. This term could come from the insertion of two twist one fields or 
one twist two field. The strong coupling contribution is again coming from the exchange 
of a single particle with a relativistic dispersion relation, last term in R2. This particle 
has mass two. We do not know whether this particle remains present at all values of the 
coupling. It might decay into two fermions as it was found in a similar context in |f22|| . 

There is an interesting interplay between the terms with integer powers of e""^ and 
the integral terms. As we analytically continue r and a, there are poles that cross the 
integration contours and we can get extra contributions which contain integer powers of 
e~'^ . In fact, they combine with the terms already present to ensure that the expansion 
has the right properties. For example, the term in R2 going like e"^"^ is not symmetric 
under a — t- —a. This lack of symmetry is cured by the lack of symmetry of the integral 
which is introduced by the iO prescription. There are similar effects that occur when we 
go to large a and then to a > r. The result is invariant under the "Wick" rotation which 
exchanges r and a. 

4.2. Expansion of the hexagon at weak coupling 

When we expand ( p.8| ) at weak coupling we generate terms going like A"r'^e~'^, k < 
n — 1, where we took into account that C starts at order A. These are the usual logs 



subtraction term. The term R^ is the contribution from the propagation of the corre- 
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that we get in perturbation theory and arise from the expansion of the exponent e(p, A) = 
1 + A71 {p) + ■ ■ ■, where 71 (p) is the one loop anomalous dimension for the excitation. More 
concretely, consider the lowest order term, which is proportional to cos (f). 



R = cos<p dp e'P^ 



(AC(i) + A2c(2) + 



■ e 



-(A7i + -)r_(_X + A2^)e-V 



(4.3) 



where we have expanded everything in powers of A and C*^*-*, are functions of p. Ti are 
the coefficients in the expansion of the cusp anomalous dimension. Now, the fact that the 
remainder function is zero at order A implies that C^^-* = c*^(p), which is the result for the 
U{1) theory, or the result at one loop. 

We can now compute the terms of order A^. Expanding ( |4.3| ) we find 



R = cos (he' 



lire' + ^c°) 
J- 1 



(4.4) 



Thus, we can take the two loop expression and look at the term Te~'^h{a). If we fourier 
transform this term we get — 7i(p)c*^(p)- We can independently compute c*^(p) from ( |3.13|) . 
We can also compute it with a trick: we can go to strong coupling, where the real result 
and the U(l) result do not mix. The U{1) result is contained in the term Ri. Its fourier 
transform is 



e'{p) oc J (icre*^'^ (cosh a log [2 cosh a] — crsinhcr) = 



TV 



p2) cosh^ 



(4.5) 



It is easy to confirm this result from the direct expansion of rf/(i-) in ( |3.11| ) at the leading 
order, see appendix D. 

We also need to know 71 (p) . This is given by the one loop anomalous dimension and 
for an F^i propagating on a sea of derivatives. It was computed in [|2^ and is given by 



(4.6) 



From |2^] (equation (3.37) with s = 3/2) we might naively expect —2-0(3) instead of 
— 2'i/'(l). However, we need to add the anomalous dimension per excitation 7^ = 3 of the 



state tr (-F+i) [|2^ to the result of [|3 
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Ho-) 




Fig. 9: Using the expressions from |21| we computed the functions h{a) appearing 
in the expansion Te^^h{a) of the two loop result for the hexagon Wilson loop. The 
dots in this figure are the values obtained by fitting the numerical data in the large 
T regime. The red curve corresponds to the analytic prediction discussed in the 
main text. Recall that there is no fit of any parameter in this comparison; hence 
this check is a very strong check of our predictions. 



We have checked numerically that this prediction is indeed true using the expressions 



from ||2^] for the two loop remainder function for the hexagon. In other words, the fourier 

in 



transform of — 7i(p)c^(p) gives a function h{a) i which is precisely the one appearing 
the term Te~'^h{a) in the expansion of the two loop result for large r, see figure fig. 9. 

It should be clear that a similar prediction can be formulated for a general polygon, 
simply by replacing C^^\p) with the leading part of logr[/(i) for that polygon. It would 
be interesting to check this prediction with an explicit two-loop calculation: the U{1) 
coefficient c'^(p) should capture the full, intricate dependence on the Wilson loop cross- 
ratios of the term re~^/i(o") in the expansion of the two loop result for large r. 

Clearly, at higher loops we can make further predictions. For example, at / loops we 
expect that the remainder function has a term of the form 



R ^ cos<Pe-^^-^f-^^ I dpe'P^c\p)Mp)f (4.7) 



This function can actually be computed analytically: 

h{a) oc cosh a (2 log(l + e^'^) log(l -|- e^^*^) — 4 log(2 cosh a)) -|- 4cr sinh a 
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This is the highest power of r that appears at this loop order. We also have terms T^e~'^ 
with k < I — 1. These also also fixed by the knowledge of the dispersion relation at higher 
loops, which is given in 0. 

5. Higher order predictions 

In principle one can push the expansion of the remainder function to higher order in 
e~'^ , and consider terms which correspond to the propagation of several particles. The 
states which appear in the U{1) answer are naively only "one particle states", cor- 
responding to insertions of a single four dimensional field on the sea of derivatives. At 
the 1-loop, in the spin chain language, such naive "one particle states" will mix with the 
continuum of true multi-particle excitations, and possibly to actual single particle bound 
states. Even without computing the exact 1-loop overlaps, we still get some constraints on 
terms of the form Te~^'^ H{a) in the expansion of the two loop result of a general polygon 
for large r, at least as long as k is sufficiently small compared to the number of sides of 
the polygon: the full functional dependence of H{a) on the cross-ratios is captured by 
the twist k coefficients Cl{p) = ^o^^pco bottom^ multiplied by some unknown function of 
p which does not depend on the specific shape of the top and bottom Wilson loops. 

6. Conclusions 

In this paper we have derived an OPE expansion for polygonal Wilson loops with light- 
like edges. The OPE expansion is performed by picking two non-consecutive null lines in 
the polygonal Wilson loop. This divides the Wilson loop into a "top" part and a "bottom" 
part with states propagating between the two. The state that propagates contains a fiux 
tube going between the two selected null lines. The states consist of excitations of this 
fiux tube. These states can also be understood as excitations around high spin operators. 
The spectrum of states is continuous and consist of many particles propagating along 
the fiux tube. In A/" = 4 Super Yang Mills these particles have a calculable dispersion 
relation . Then the OPE expansion leads to predictions for the Wilson loop expectation 
values. Namely, it implies constraints on the subleading terms in the coUinear limit of the 
remainder function, which is the function containing the conformal invariant information 
of the Wilson loop expectation value. We have checked these predictions both at strong 
coupling and at two loops at weak coupling. 
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In order to perform this expansion at strong coupling, we have derived an alternative 



presentation of the integral equations p5|,|T^,0 . In the new representation the parameters 
in the equations are the physical spacetime cross ratios. In addition, the formula for the 
area is given by the Yang- Yang functional associated to the modified TBA equations. 
This new form of the equations are as useful as the old ones ||2^,|T5|,|T^ for numerical 
computations. 

For single particle exchanges we have been able to characterize completely the state 
of the particle in terms of its momentum. Using the exact formula for the single particle 
dispersion relation, given in one can find an exact prediction for the first subleading 
term in the coUinear expansion. Expanding this formula in powers of A one would obtain 
specific predictions to all orders. Here we have checked this prediction with the only 
available direct computation in the literature, the two loop hexagon [^^Jl^,^,^ . 

When we have multiparticle exchanges the situation is less clear and we will prob- 
ably need the full power of integrability, together with the infinite number of charges to 
characterize the state. It is likely that characterizing the states in this fashion and then 
demanding a consistent expansion in all channels one could determine the full function. 
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Appendix A. Describing the family of polygons in terms of momentum twistors 



One can describe the kinematics of the family of polygons that we considered in 
section 3.1 by twistors. These are sometimes called "momentum twistors" and they appear 



naturally both in weak coupling [26| and strong coupling [25.15,16| calculations. 
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Each vertex of the polygon can be represented by a null vector Z in i?^'^, defined up to 
rescaling. The vector Z can be thought of carrying two antisymmetric spinor indices, and 
can be rewritten as the bilinear antisymmetric combination of two twistors. Schematically, 
Z = A A A'. Both A and A' carry a spinor index of i?^'^ and are defined up to rescaling. 
A polygon with N null sides can be given as a sequence of N twistors A^, such that the 
intersection of the sides % and z+1 is the point = AjAAj+i. Indeed Zi_i/2-Zij^ii2 = 

is the condition that the i-th. side should be null. It is also useful to introduce dual 
momentum twistors fXi, cospinors which can be defined as fXi = Aj_i A Aj A A^+i, and also 
safisfy Zi+i/2 = fJ^iA fJ^i+i- 

A twistor A together with an orthogonal dual twistor (//, A) = 0, defines a null 
line (and viceversa): all points of the form Z = X Av with = 0, v defined up to 

rescaling and shifts by A. Our family of polygons with fixed sides i,j must have fixed 
Xij and Hi J. We are now in condition to give an alternative description of the group of 
transformations fixing two non-intersecting null lines, defined by the pairs Aj, Hi and Aj, 
Hj. In order for the lines to be at generic position with respect to each other, {iJ,i,Xj) 
and (/Xj, Xi) should be both non-zero. A conformal transformation M will leave the lines 
invariant if Xij are right eigenvectors of M, i.e. MXi = e^'^'^Xi and MXj = e~"''^'^Xj, and 
Hi J are left eigenvectors. As long as {nijXj) and {iJiji\) ^tre non-zero, it must be that 
HiM = Hte~"^^ and HjM = Hj^"^^ ■ 

If we decompose the twistor space into multiples of A^, multiples of Aj and the or- 
thogonal to M will be block diagonal, with elements e^^^'^ and e~^R for some SL(2) 
transformation R. To pick a reference square is the same as picking two more twistors A± 
in the orthogonal to Hij- The square is left invariant by M iff X± are the two remaining 
eigenvectors of M, with eigenvalues e^'^~^ . 

Once we have picked a specific reference square, the family of polygons is defined by 
the twistors A^ • ■ • Xj and MXj-^i ■ ■ ■ MXi-i. 



Appendix B. The hamiltonian picture and its analytic continuation 

In this appendix we consider the large spin limit of local operators and we see how 
they give rise to the flux vacuum and its excitations. 

We write the coordinates of R^'^ as Zm = Z_i, Zq, Zi, ■ • ■ Z4, with the condition 
= 0, where the indices are contracted with the Minkowski metric of R^'^. In addition 
we impose the identification Z ~ XZ. Usual poincare coordinates are 

a;^ = — , for /x = 0,l,2,3 (B.l) 
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The usual R^'^ metric is simply the induced metric on the lightcone = with the gauge 
condition Z_i + Z4 = 1. Different "gauge fixing" conditions lead to Weyl transformations 
for the metric. Just for reference, the usual coordinates of R x are defined by 

Z Z^ 
tant = — ^, n' = ^ = , z = l,---,4 (B.2) 

where n* is a unit vector in i?^ and describes a point on . Similarly, if we choose a gauge 
fixing function Z| + Zf = 1 we get a metric which is that of AdSs x 5"^. This is again 
conformally related to R^'^. 

Within this AdS^ factor we can choose coordinates which are similar to Euler angles 
for (since AdSs is an analytic continuation of 5"^). We can write 



'2 "T ^3 

where are the usual Pauli matrices. We have parametrized the space in terms of coor- 
dinates Ti, (3, a and (j) = arctan(Z'2/Z3). The metric in these coordinates is 

dsAdSsxS^ = - rfrf + rfcr^ + + 2 sinh 2l3dTid(j + dcj? = 

^{da + sinh 2l3dTif - cosh^ 2l3dTf + d/?^ + dcf)"^ ^^'^^ 

These coordinates make manifest the symmetries of the two Wilson line configuration. 
We see that the cr direction is fibered over an AdS2 space which realizes the 5'L(2, R) 
symmetry. This metric differs from the fiat metric of R^'^ by an overall conformal factor, 
which is not important if we are dealing with a conformal field theory. We can consider 
two Wilson lines at Z2 = Z3 = and at Zq/Z^i = Z1/Z4. The SL{2^R) that preserves 
two of the lines acts by multiplication on the left in ( p.3| ). These Wilson lines lie at 
the boundary of the AdS^ space and produce a flux of color electric field along the Ti,a 
directions, F^-^fj =constant. This is a flux tube whose energy is localized in the transverse 
non-compact (3 direction. This is discussed in more detail in [jT^. Excitations of this flux 
vacuum constitute the states that arise in the OPE expansion. In any theory where the 
flux cannot "break" , these states are well defined. For example, in A/" = 4 super Yang Mills 



This AdSs has no relation to the AdS space that appear in AdS/CFT. This AdS space is 
purely within the boundary theory. This discussion applies to any CFT, whether is has a known 
gravity dual or not. 
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this flux is well defined at finite for any value of the coupling if the external lines are 
in the fundamental, so that the flux cannot be screened. In theories with fundamentals, 
these states might be not deflned away from the planar limit. 

Let us now discuss the analytic continuation r; — )► zr. We also need to set /3 — )■ 
Then the coordinates ( p.3|) continue to describe an AdSs space written as an SU{1,1) 
group element where we have analytically continued — )■ iYo and Yq — ?■ zl^. This now 
maps the Wilson lines as follows. After the analytic continuation the parametrization ( [B.3| ) 
becomes a parametrization as an SU{1, 1) matrix. We can perform a simple relabeling of 
the matrices a2 — > o"3 and as — ?■ — a2 which corresponds to the transformation which takes 
the SU{1, 1) group element back to the usual presentation in terms of SL{2,R) matrices. 
After this transformation we end up with a parametrization of the form 

We now see that the transformation r — > r+constant corresponds to a combination of 
an ordinary dilatation and a boost if we go back to the usual Poincare coordinates via 
( p.l| ). This is such that it leaves a;"*" flxed and it changes x~ . We then get a square as in 
flg. 4(c). The two original Wilson lines are the ones acted on by these shifts of r and are 
the ones that are extended along x~ . The other two sides of the square have arisen after 
the analytic continuation and they are the lines on which are are inserting the operators 
creating the flux vacuum together with its excitations. If we only have the flux vacuum, 
these extra lines do not have any further operators inserted. See also flg. 6(b). 

B.l. Analytic continuation at strong coupling 




Fig. 10: Plot of the the spinning string solution for small spin and for large spin. 
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Fig. 11: Plot of the Euclidean continuation of the spinning string solution for 
small spin. Side view and top view. 




Fig. 12: Plot of the Euclidean continuation of the spinning string solution for very 
large spin. The edges of the string in the bulk approach the boundary, leading to 
the two other boundaries of the Wilson loop. 



At strong coupling we can consider the string solution that is dual to the high spin 
operators considered in 0. In this appendix we will show how the high spin limit, plus an 
analytic continuation, produces the square Wilson loop. A closely related discussion can 
be found in . The idea that the large spin limit can be described in terms of Wilson 



loops can found in ||2^. It is convenient to focus just on an R^'^ subspace of the boundary 
theory and an AdSs subspace of the bulk, with metric ds'^ = — cosh^ pdt^ +smh.^ pd(f'^+dp'^ 
This should not be confused with the AdSs subspace discussed above which was purely on 
the boundary. For finite spin we have the solution 

(p = wt, p e [0,pmax] , tanhp^ax = l/w (B.6) 

where w is very large for small spin and w — )■ 1 when the spin goes to infinity. This is the 
Lorentzian time solution. We can now analytically continue t — > zr and f — )■ ix- We then 
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write 

y_i ±1^4 = coshp , Yo±Yi = e^'"" sinhp , p G [0, pma.] (B.7) 

This is plotted in fig. 11. We see that as we approach the boundary and r — )■ oo the 
solution approaches two null lines. The worldsheet tips, sitting at p = Pmax^ move through 
the bulk and approach the boundary along these null lines. As the spin goes to infinity 
— 7- 1 and pmax — ^ oo. Then the tips of the string get closer to the boundary. In the 
limit, the tip is joining a point that is a quarter way from the tip of the null line in fig. 11 
to a light-like separated point on the other line. 

At weak coupling we expect a similar picture. In fact, a closely related picture was 
discussed in [^. The insertion of an operator Tr[^df^] produces a displacement of 
the field insertions along the x"*" direction together with an adjoint Wilson line. As 
the spin gets larger we expect that the effective displacement of the two fields along the 
null direction becomes larger. Now suppose that we have two such operators inserted at 
antipodal points in the cylinder. These points are spacelike separated. However, as we 
increase the spin, each of the points is splitting into two and they become displaced along 
the null direction. As this displacement grows one of the $ insertions of one operator 
can become light-like separated with the $ insertion of the other. This happens as the 
spin goes to infinity. Operators that diagonalize the dilation do not have fixed separation 
along the direction, but are a suitably weighted superposition. At strong coupling this 
manifests itself as the fact that we always get full null lines as we approach r — )■ ±00 in 
O, see fig. 11. 

This picture is not restricted to the large spin limit of twist two operators, but it 
would also hold for higher twist operators which contain other field insertions among the 
derivatives, as long as we consider the large spin limit and look at the lowest excitations. 
See [0 for further discussion and further references. 

Finally, notice that by inserting n high spin twist two operators at spacelike separated 
points and by taking their large spin limit we expect to reproduce a null polygonal Wilson 
loop with 2n sides. In particular, the three point function of high spin operators of the 
type Tr[$9'^$] should produce the hexagonal Wilson loop, when each of the spins goes to 
infinity. The cross ratios of the Wilson loop should come from the orientation of the spins 
and their ratios as they go to infinity. It should be interesting to see if this produces a 
simpler way to compute Wilson loop expectation values. 

Alternatively, a null Wilson line between two $ insertions can be expanded in derivatives. 
Then, correlators of these Wilson operators are dominated by large spin. 
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Appendix C. Remarks on the dispersion relation 

In this appendix we would like to make a couple of remarks on the dispersion relation 
for excitations around the flux vacuum described in the above appendix. Let us discuss 
it flrst at weak coupling. In the free theory, the particles do not feel the flux and the 
spectrum is the same as in the absence of flux. It turns out that the energies are quantized 
and independent of the momentum, due to the SL(2, R) xSL{2, R) representation theory of 
the problem . At flrst order in the coupling we get a correction. The states in question 
can be viewed as excitations around an inflnite "sea of derivatives" . In other words, high 
spin operators have the rough form Tr[$9^$]. In the large S limit, the two field insertions 
$ become displaced by the derivatives giving rise to the Wilson lines. Thus, we are left 
with the sea of derivatives. We can then insert other fields among this sea of derivatives. In 
the free theory the total twist, which we are calling "energy", is just the twist of the extra 
fields, which could be 1,2,3, etc. At first order in the coupling we see that the degeneracy is 
broken and we get a non-trivial dispersion relation. A similar phenomenon occurs around 
other large charge "vacuua" such as the BMN vacuum ||2^. At one loop the twist one fields 
give rise to well defined excitations with a dispersion relation e{p) = 1 + X'y{p)- For A/" = 4 
SYM these excitations are two gauge field insertions F^i and six scalars and eight 
fermionic excitations i^.^.! ^. Since these fields are well defined particles at weak coupling, 
we expect that they will survive at all values of the coupling. At twist two we have several 
possibilities, one of them is for example. We do not know if it survives as a well 

defined particle for all values of the coupling. At strong coupling we have an accidental 
relativistic symmetry, as we have around the BMN vacuum [^|. This is only present at 
strong coupling and is broken as we go away from the strong coupling limit. We have ^0 



particles with mass m = \/2 that corresponds to the particles. We also have eight 
massive fermions of mass m = 1. We have five massless scalars from the five sphere. Due 
to strong IR effects these give rise to six massive particles of very tiny mass m oc e~^^ 
Finally, there is a particle of mass m = 2. 

It turns out that some of these particles can be interpreted as "massive goldstone" 
particles. 

Notice that after the analytic continuation in (|B.5| ) the coordinates ( [B.4| ) become 

rfs2 _ ^^2 + ^^2^2 sin 2pdrda - d/S"^ + d(p'^ (C.l) 

We see that we have a symmetry under t ^ a and (3 — t- — /3. This leaves the fiux 
essentially invariant. This implies that the dispersion relation is expected to have the 
following symmetry. If we write the dispersion relation as /(e,p. A) = 0, then if e and p 
are a solution, then so is e' = ip ,p' = ie. Thus we have a Wick rotation symmetry. 
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C.l. Massive Goldstone particles 

The flux vacuum that we are considering in this paper is breaking some symmetries 
spontaneously. When we have spontaneously broken symmetries we can act with the bro- 
ken symmetry generators and generate some particular excitations. These excitations have 
particular energies and momenta which are fixed by the commutation relations between 
the broken symmetry generator and the energy and momentum generators. Let us see 
what this implies for our case For example, the flux breaks supersymmetry. Under 
the SL{2)l X SL(2)ji C 5*0(2, 4) generators the supercharges transform as (^,0) + (0, ^). 
Thus they have either e = ±1 and p = Oore = 0, p = ±i (recall that the energy is 
a generator inside SL{2)l and the momentum is a generator inside SL{2)r). Thus the 
dispersion relation of the fermions should obey 

e(A,p = 0) = l , e(A,p = ±z)=0, fermions (C.2) 

This should be true for all values of the coupling. E.g., this argument determines the mass 
of the fermions at strong coupling. 

For the bosons associated to Fj^i there is a similar prediction. In this case the corre- 
sponding generators have spin (i, i) under the SL{2) x SL{2) of AdS2,. Thus, in this case 
the Goldstone bosons have e = ±1, p = ±z. Thus, in this case we have the constraint 

e(A,p = ±z) = 1 , bosons (C.3) 

One is tempted to make the same argument for the broken generator involving the broken 
SL{2, R) generators. In this case one would be lead to a particle with energy e = 2, p = 0. 
This is indeed present at strong coupling, it is the m = 2 particle. On the other hand, 
it might be that this particle decays into pairs of fermions, which would mean that the 
symmetry generator would create pairs of fermions rather than a single boson. Notice that 
the prediction ( |C3| ) is quite general. It should hold in any CFT in dimensions d > 3. 
Superficially this relation seems to be contradicted by ( [4.6| ); we believe that the limit 
p ^ i does not commute with perturbation theory. A closer analysis of the finite coupling 
dispersion of should hopefully confirm this. 

This non-relativistic version of the Goldstone theorem constrains the dispersion re- 
lation only at particular values of the momentum. In these cases the broken symmetries 
do not commute with the Hamiltonian and thus they give rise to excitations with finite 
energies. 

Notice that the excitations around the flux vacuum have a gap, both at weak and 
strong coupling, though the gap becomes non-perturbatively small at strong coupling. In 
particular, the anomalous dimensions of the $^ insertions at zero momentum is negative, 
as naively expected from the fact that they become very small at strong coupling p^,p|. 
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Appendix D. On the expansion for Wilson loops in the U{1) theory 

In this appendix we aim to give some more details about the computation of ^[/(i) for 
a general polygon. It is convenient to rewrite the contour integrals in terms of projective 
coordinates Z as 

log.„,, ^ / / d.' (D.l, 

^top ^bottom 

Notice that a position-dependent real rescaling Z{s) — > Z{s)X{s) changes the integrand up 
by a total derivative only, which integrates to zero. Be A,D the top vertices of the square, 
B,C the bottom vertices. Wsquare — Wtop consists of a path 7 going from Ato D minus the 
straight segment AD, see fig. 8(b). Wsquare — Wbottom consists of a path 7' going from B to 
C minus the straight segment BC. We can write 7(5) = a{s)A + b{s)B + c{s)C + d{s)D + 
n{s), where a,b,c,d are scalar functions and n is a vector in the plane orthogonal to the 
square. We can also write Y{s') = a'{s')A + b'{s')B + c'{s')C + d'{s')D + n'{s'). It is useful 
to parameterize the segment AD as Zo{s) = a{s)A + d{s)D and write 7(5) = Zo{s) + 5{s). 
Without loss of generality, we can pick a{s) = and d{s) = e^. We will also parameterize 
the segment BC as Z^(s') = b'{s')B + c'{s')C and write Y{s') = Z'q{s') + 5'{s'). Without 
loss of generality, we can pick b'{s') = e~*' and c'(s') = e*'. We can also normalize 
A C = B ■ D = 1/2. 

Notice that 5 ■ Z' = 5' ■ Z = and their derivatives are also orthogonal. After some 
simple algebra 
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The action of M on 7'(s') can be combined with a trivial rescaling and shift of the 
integration variable as e-''MY{s' -a) = e'^'^-'^a'is' -a)A + e-''' B + C + e-^''+'^d' {s' - 
cr)D + e~'^n'^{s' — a). We denoted the action of the rotation in the transverse plane on 
n' as n'^. Hence in order to write the OPE we want to redefine 5'{s' — a,T,a,(j)) = 
e-^^-^'a'is' - a)A + e'^^+^d^s' - a)D + e-^n'As' - a), and 
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, , , ,5{s)-5'{s' -a) + 5{s)-5'{s' -a) , 

logr[/(i) = / ds ds' ^ \ , „ ' \ \ ^ D.3 

^ ' I I cos\v{s - s') + 5{s) ■ 5' {s' - a) ^ ' 



—00 •/ —00 



32 



Now we can expand the denominator and Fourier transform term by term. The leading 
piece is 



^n{s) ■ n'^{s' -a) + h{s) ■ n'^{s' - a) 



/oo poo 



+ 1 



For the hexagon we find that n{s) ~ e~l^lno where uq is a constant vector. We have a 
similar expression for n' . Its Fourier transform is then proportional to (1 We get 

a similar factor from n' . Together these to factors give the function in (^1.5|). 



Appendix E. The modified TBA and Yang- Yang functional for AdSs null poly- 
gons 

In this appendix we present a new version for the integral equations that determine 
the result at strong coupling in the case that the polygon can be embedded in R^'^ and 
the surface in AdS^. This modified integral equations, or TBA equations, involve only the 



physical spacetime cross ratios, in contrast with the ones in [16| which involved some other 



auxiliary parameters. In addition, we will find that the expression for the area can be 



written as the critical value of the associated Yang- Yang functional [pl| , |3^ . The compu- 
tation of the regularized area (or better, the remainder function) for the minimal surface 
in AdSs ending on a given null polygon on the boundary is done in |2^,|l6| through three 
basic steps. The first step is to promote the cross-ratios of the null polygon to holomor- 
phic functions of a spectral parameter which capture the higher conserved charges of 
the classical integrable system. The second step is to derive a set of functional equations, 
which are then converted into convenient integral equations. The third step is to compute 
the regularized area from the higher conserved charges or ( dependent cross ratios. 

The specific problem of minimal surfaces ending on null polygons in AdSs is a special 
case of a general theory of Hitchin systems on a Riemann surface. The results of this 
appendix apply to any such system, hence we will try to keep our integral equations as 
general as possible. Fortunately, the integral equations given for general Hitchin systems 
3^,0 and the integral equations which are optimized for the case of polygons in AdSs 



llql actually coincide in the simplest kinematic region. There is a large amount of freedom 
in setting up the Riemann Hilbert problem, which leads to a variety of different forms for 
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the integral equations, each with its own advantages and disadvantages. The purpose of 
this appendix is to present the integral equations in a way which is well suited for the 
analysis of soft limits, and make manifest the cancellation of spurious terms between the 
various contributions to the regularized area. 

The TBA-like integral equations for general Hitchin systems can be written as 

lnX,(C) = |^ + z^, + Z,C-T^$^0(y)(7,7) / f^^log(l + Xy(0) (E.l) 

The equations depend on some discrete data: the set F of possible labels 7, 7', certain 
integer numbers ^(7), an antisymmetric pairing (7, 7'). This data is the final product of a 
careful WKB analysis of some differential equations on the surface, and for the purpose of 
this appendix we only need to use some basic facts about it. The are auxiliary complex 
numbers which we wish to eliminate from the equations. The are angles which will 
presently be set to zero (or tt, in which case they can be reabsorbed by changing the sign 
of some X-y in the equations): in the language of Hitchin systems, we are restricting to a 
real section. The lines of integration are straight rays from ^' = to C = 00, which for 
convergence reasons should lie in the half plane where Z^/C,' has a negative real part. A 
canonical choice is to set them to Z^/C' real and negative. The relative ordering of the lines 
is important, in the sense that if two lines are moved across each other, the integration 
contours will cross poles of the integration kernels. The result of such a "wallcrossing" is 
captured by a certain specific change in the discrete data. 

Once the angles 6*-^ are set to zero, a simplification occurs: C) = X^{C,). This 

is due to a certain symmetry of the equations: labels in F come in pairs, which we can 
denote as 7 and —7, such that 0(7) = 0(— 7), (—7,7') = (7, —7') = — (7,7')5 -^-7 = —Z^. 

We can combine the contributions from 7' and —7', and have the sum run only half 
of the set of 7's, say F+. From now all sums will be over this subset. We will also write 

\r.X,{e) = Z,e-' + Z,e' + — 0(7')(7,7') / • Jn, ,^ og{l + X,,{e')) (E.2) 

In order to make contact with the equations of for 2N gluons we can take the 
labels 7 = s to run over integers 1 ■ ■ ■ — 3, with = 1 and (s, s + 1) = —1. We also 
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need to set Zg = —^rris, and use the canonical choice of hnes. The functions Xg and Yg 
agree on the integration hnes, up to an appropriate shift of 9 by argZ^-Ill 

We want to rewrite the equations in terms of and x~ , the values of the cross-ratios 



Xj at the physical values C = 1 and C = ^- This is easy: we set = 1 and ( = i in ( [E.2|) , 
solve for Zj and Z^ and then insert the result back into ( [E.2| ). We find 

In X^{9) = cosh 9\nx^ — i sinh ^ In + 

1 r^. /^/ /\ /" d^' sinh 2^ , ,^ ^ (E.3) 

Notice that the auxiliary parameters Z^ have disappeared from the equations and we 
only have the spacetime cross ratios x^ . Notice that the kernel of these modified TBA 
equations is not symmetric in 9 and 9'. It becomes symmetric if we pick a different choice 
of rapidity variable, u = gf^h 29 ' that —2d9' = du' sinh^ 29' . The usefulness of this 
alternative rapidity variable will become evident momentarily. General TBA equations 
with a symmetric kernel 

lnXa{u) = La{u) + ^y] [ du'Kab{u,u')log{l+Xi,{u')) (E.4) 
can be recast as the conditions for a Yang- Yang functional to be extremized 

+ ^y^f du I du'Kabiu,u')pa{u)pb{u') 

8^ ^ Jia Je. 



(E.5) 



Indeed the variation with respect to p sets 

<^a(w) + 77- / du'Kabiu, u)ph{u) = (E.6) 
271- V J I, 



b -""b 

while the variation with respect to (pa sets 



Pa{u) = log + e^-H-'^"(-)) (E.7) 



The reader should be cautious in extending this relation away from integration lines: the Ys 
functions are usually defined as analytic continuations from the integration lines, while the Xs 
are defined by the integral equations for all 6 and therefore have discontinuities. 
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If we set log Xa = La— (pa we recover the TB A equations (|E.4| ) . We would like to show 
now that the interesting part of the regularized area (or better, reminder function) of the 
minimal surface in AdS^ coincides with the extremum of the Yang- Yang functional for the 
modified TBA equations. In a sense, this result makes manifest an important property of 
the area: it should be the extremum of an action functional with fixed boundary conditions 
given by the choice of physical cross-ratios. 

For convenience, let's take to be odd, where the polygon has 2N sides. For even 
N there are some slight complications, which pop out in various places in the calculation, 
only to cancel out at the very end in the remainder function. For that reason, it is simpler 
to just do formal computations for odd A, and then recover the even A results by a 
soft coUinear limit. If A is odd, there is an antisymmetric matrix w-y,^' which roughly 
speaking inverts (7,7').lll More precisely, it satisfies 7' -^7^7,7' (t'' 7") — ^1" 

Ey,7"(7,7")«'7',7"(7",y") = (7,y")- 

The higher conserved charges of the classical integrable system are hidden in the large 
positive 9 asymptotic expansion of log A^ ~ Xl^-i ^71,76""^. An alternative set of charges 
appear at large negative 9: log A-^ ~ J2'^=-i £71,76"^^. The contribution to the regularized 
area denoted as Aperiods + A free in can be computed from the conserved charges as 

i y w^^-y'C-i^^ci^Y . An alternative expression is —i J^-y 7' w^^yC-i^jCi^-y' . The cleanest 
formulae usually arise by averaging the two expressions. 

From the TBA equations (|E.3|) we can compute 

c_i,7 = i(lnxj-^lnx-)-— 5^ 0(y)(7,7') / " h2^ + (E.8) 



1 /■ r]9'p'^^ 

ci,7 = i(lnx^ + zlnx,-)-— 0(y)(7,7') / + (E.9) 

c_i,7 = i(lnx++^lnx-)-— 5^ 0(7')(7,7') / + (E.IO) 

7'er+ ''^1' 

ci,7 = i(ln<-zlna:;)-— j ^ ^^^^""^^^ + ^^'^^'^^ (^.11) 



7'er+ 
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In some cases the set of 7's is an over complete basis, that is why we cannot find a proper 



inverse. In the kinematic region considered in [16| one can find a proper inverse. 

Recall that A denotes difi^erent pieces of the area. The Wilson loop expectation values are 



R 



then obtained by {W) ~ e"^^^'''^"^ ~ e-^^^""""). 
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The piece of Ap^riods + A free with no integrals is 

Aq = -^'^Wj^ylnx^lnx~, (E.12) 



7:7 

The piece with one integral, averaged, is 



Atenip = -^ 5Z / ^siTh2^^ (sinh6'lnx+ - z cosh 6* In ) log (1 + X-^(6')) 

(E.13) 

Notice that the term in parenthesis is dg (cosh 9 In — i sinh ^ In x") . We can use the mod- 
ified TBA ( |E.3|) to trade it for de InX-y(^) up to terms which we will combine with the other 
two-integral pieces of Aperiods + Aj^ee- Then we can write dglnX^iO) log(l -\- X^i{6')) ~ 
de {—Li2{—X^)) and integrate by parts to 



Ai = -y: fi(7) / -^Lz,{-X,) (E.14) 

smh 26 



■yer+ 

The remaining term is 

de' r de 

7,7'er+ 



sinh 26* 



X cosh20a, ^-^^^J log (1 -f (^')) log (1 + x^ie)) 



(E.15) 



We can symmetrize the kernel to sinh(e'-e) +sinh(^^— ^) cosh 2 (^-(-^'). The second piece 
cancels out against the remaining two integral terms and we are left with Aperiods+Af^ee ~ 
AQ + A1 + A2, where 

.inh2».inh(e'-9) 

log(l + Jfy(e'))log(l-fX,(9)) 

We see that Aperiods + A free = ^0 + ^^cr, where YYcr is the critical value of the 
YY functional. In summary, the full area of the surface can be written as A = 
Adrv + ABDS-Hke + Aq + YYer, where YYer = Ai + A2 m (FtI and ( WM - are 
the divergent terms and AsDS-uke can be found in formula (5.10) in [^. We see that 
YYcr is manifestly small in the region where the X^ are small, which is when the x^ are 
small. We also expect that AsDS-iike — Abds + Aq becomes small in this region. This is 
an expression written purely in terms of the spacetime cross ratios. 
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(E.17) 



E.l. Evaluating the OPE for general configurations in AdSs 

In this subsection we expand the OPE with AdS^ kinematics in the case that we spht 
the polygon into two general subpolygons. 

If we are given any functional F[x] and a small deformation F[x] + f[x,y], we can 
ask what is the relation between the critical values Fcr and {F + f)cr- It is elementary to 
show that the difference of the two critical values can be computed as the critical value 
of /[xcr, y], where Xcr extremizes F{x). We can immediately apply this to the case of the 
Yang- Yang functionals. We take F = YY[c\, the Yang- Yang functional with a certain label 
c G r+ erased, and / is the remaining part (notice that kernel Kcc is always zero for us) 

/" rfw(pe(^l)(/>e(^/)-L^2(-e^^(")-'^^^")))^- 

+ -r-^y^ I du I du'Kcb{u,u')pc{u)ph{u'] 

The extremum of / with respect to pc and (pc is simply 

U = YY^r - YY[c],r I duLi2{-X,{u)) ~ 7^ /" duX,{u) (E.18) 

where 

lnX,{u) = L,{u) + ^y^ [ du'K^hiu, u') log (1 + Xtiu')) (E.19) 

The soft limit OPE we consider in this paper leads exactly to this sort of decoupling 
limits. For simplicity, we can consider the simple kinematic region of In the AdS^ 
case, if we decompose a 2A^-gon into a 2?! + 2-gon and 2N — 2n + 2-gon by applying the 
r rescaling to the x+ coordinates, it is easy to see that Z^-i acquires a large real part, 
so that x^_i becomes very small, scaling as e""^, while all the other x^ and x~ are held 
fixed. Indeed, — 1] is the sum of the Yang- Yang functionals for the 2n + 2-gon and 

2N — 2n + 2-gon. Then the leading exponential behavior is controlled by 



U = YY,r - YY[n - / rfwe'"""-i ^^^^^--^^^^-i Vi(6')C2(6') (E.20) 

where 

\nCi{e) = ^ V / du'K,,{u,u')\og{l + X,{u')) (E.21) 

and 

lnC2(^) = ;^ V / du'K,,{u,u')\og{l + X^,{u')) (E.22) 

s>n-l ^'^b 

are interpreted naturally as the density of one-particle excitations on the GKP string 
created by the 2n + 2-gon and 2N — 2n + 2-gon respectively. 
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Appendix F. The modified TBA and Yang- Yang functional for AdS^ null poly- 
gons 

In this section we generalize the results of the previous section to the most general 
AdS^ kinematics. For convenience, we take the number of gluons n to be odd. The case 



were n is even can be recover by a soft coUinear limit. In the notations of [jT^ the TBA 
equations take the form 

log I2, 3(6*) = - 1 I a/2 cosh (6' - i(j)s) - K2 as - Ki -k I3s 

logYi^siO) = -\ms\ cosh{e -i(j)s) - Cs - \K2-k Ps - i^i ★ - (F.l) 
logY3,s(6') = -|ms| cosh(6' -i(t)s) + Cs- \K2-k Ps - i^Ti ★ + 

where -k denotes convolution. The kernels are 

K^{e) = , K2{9) = ^^^^ , Km = Atanh2^ (F.2) 

27r cosh i) TT cosh 26' tt 



and 



as . log , a + ^M) (; + >-3.) ^ (1 + Y.^s-.) (1 + 



a + Y2,s-i) (1 + Y2,s+i) ' ' " (1 + (1 + Y^,s-i) 

°^ (1 + (1 + (1 + n,.-!) (1 + i"3,.+i) ■ 



(F.3) 



We now have three parameters per column s: the magnitude of the mass |ms|, its phase 
(^s and the chemical 'potential Cs- As explained in the AdS^ case we can eliminate them in 
favor of physical cross-ratios ya,s = Ya^s{^) where a = 1, 2,3. The hatted F-functions are 



defined in the appendix D of IT^] and are given by Ya,s{0) = Ya^siO) if a + s is even and 
Ya,s{d) = Ya^s{d — ^7r/4) for + s odd. The modified TBA equations then read 

log 1^2, s -Es = -K2 o&s - Kio Ps , 

log Yi,sYs,s - V2Ei^-^'^'^ = -Kf (-)'^'l o 4 - 2Ki o , (F.4) 
^ogYi^s/Y3,s - log 2/1,5/2/3,5 = -^3 7s , 

where f^'^^{6, 9') = f{6+imn/4, 9'), {dg, /3s, 7s} are {os, /^s, 7s} evaluated on the hated F's 
and Eg = (v^sinh[6' + (—)^^7^/4] log 2/2,5 — sinh[^] log yi,s2/3,s) is the new source 

which depends only on the physical cross ratios. The modified kernels read 

jy _ 1 sinh(2g) y _ i sinh(2g) 

^ 27r sinh(2^') cosh(^ - ^ tt sinh(2^ - 2^') smh(2^') ^ ' ' 
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and 

K2 = zy2sinh[^ -e' + {-yin/4]Ks . (F.6) 

Notice that in ( [F.4|) , o no longer denotes a convolution but instead an application of an 
integral kernel K{6,6'). Let us denote the right hand side of the three equations in (|1^\4| ) 
by ^2,s, Ai^s + ^3,s and Ai^s - ^3,s respectively. The contribution Aperiods + A^xtra can 
then be computed as explained in the previous sections. The result turns out to be simply 
given by 

■^periods ~l~ -^extra ^0 YYq , (F.7) 



where 



de 



7rsinh^(2^) 



Li2{-Ya,s) - \ log(l + Ya,s)Aa, 



(F.8) 



is the value of the Yang- Yang functional at the extremum. Note that the natural particle 
rapidities are again given by u. 

The pice Aq is very much like Ap^riods and is given by Aq = -^Va,si^a,s-a's'Va's' where 



= - [2 logy2,s - (1 + ti-ir+') logyi,.2/3,.] (-1)^+^1 - {-ly+'tY-^ , (F.9) 



Va,s is the complex conjugate vector, and uja,s;a's' is the inverse of the intersection form of 



cycles and is given in [16 



F.l. Explicit details for the hexagon 

In order to illustrate how the kernels of the modified TBA and the Yang- Yang func- 
tional emerge in a concrete, relatively simple example, we can consider the case of the 
hexagon in full, tedious detail. We are interested in a limit where one of the three cross- 
ratios, conventionally U2, goes to zero, while the other two remain finite, and their sum 
Ui + Us approaches 1. At strong coupling, this conditions constrain the behavior of the 
conserved spin four current P{z): the two zeroes of the polynomial P{z) are far from 
each other in a specific direction. In the language of [jl5| this limit corresponds to a large 
absolute value of the "period" Z with phase ip ~ — 7r/4, (^ = (/?-|-7r/4~0 (more precisely, 
we will take it to be slightly positive). 

Compared to the previous section, we now have only one value of s and only three Y 
functions Ya = Yi^a- One of the three TBA equations is trivially solved, since the ratio 
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Yi/Ys is constant. The other two equations for logl2 and logYiY^ reduce to the equations 
in |T5| (see (3.6) and (3.7) in |T|]) 



£(6* - z^) = ze'"/^-^ + ze^-^"/4 + j de'K2{e-e')L{e') + J de'Ki{e-e')L{e') (f.io) 

e{e-iifi) = V2(Ze*"/4-^+Ze^-*"/4)+ J de'2Ki{9-e')L{e')+ j dO' K2{e-e')L{e') (F.ll) 

We denote as L(^) = log(l + e-^'^^"*"^)) and L(^) = log(l + /xe-^(^-^'^))(l + ^u-^e-^^"^'^)). 
We shifted the integration variables for convenience. Remember the correct integration 
paths: the real part of 6' runs from — oo to oo, the imaginary part is fixed at (p. 

We can evaluate the equations at = and in /A in order to trade Z for cross-ratios. 

Z + Z = log(6i) - j de'K2{e' - m/A)L{e') - j de'Ki{6' - m/A)L{e') (F.12) 

and 

v/2(Ze*"/4 + Ze"^"/^) = log(l/w2 - ^) - j de'2Ki{e')L{e') - j dO' K2{9')L{e') (F.13) 

We need to decompose 
^gZ7./4-e^£g9-i^/4 _ _^(z+Z)y2sinh(^)+^(Ze*^/^ + Ze-^'^/^)y2sinh(^-i7r/4) (F.14) 
so that we are ready to substitute 

ze-/'-' + Ze'-I^ = E{e) [ ^g^^^^M^ + ^0 L(^0 - 1 / de''""'^^^' ^^{9') 

^ ^ n J sinh2^' ^ ^ n J cosh 2^' ^ ^ 

(F.15) 

in the TBA equations. Here we define 

E{9) = -zV2sinh^log(6i) + z sinh(0 - z7r/4) log(l/u2 - !)• (F.16) 
We get the new kernels 



^/.x V2 f cosh2^sinh(^-^') ~ 

-J- fd9' r(9') 

2n J cosh(6'-6'')cosh26'' ^ ^ 



(F.17) 
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and 

6(6 - 10) = V2E{e)-- f dd' ^-^^l^-——L{e') 

_V2 f sinh2^sinh(g-g0 ^^'^^^ 
TT J cosh2(6'-6'')cosh26'' 

In order to compute Aperiods + A free it is useful to go back to the original expression 
for the regularized area, involving the conserved quantities hidden in the large 9 asymptotic 
expansion of the e and e functions. If we denote the coefficients of the expansion as 

oo oo 

6(0) ~ Yl ~ E (^-19) 

n=—l n=— 1 

then 

Aperiods + Afree = ~ {e'^^^e^ih - c'^/^e^^e^^ . (F.20) 
We can expand in powers of 



e 



361' 



e(e - i<fi) ~ E(e)-— / d9' -(^^ + — — )e^ e'^ 1(9') 

^ ^' ^ ' J l^4^sinh^' cosh^'^ 2sinh2^' ' ^ ' 

1 f e'+'' e^^' 



(F.21) 



"2^ j [^^^ " ^^^^2^'" j ^^^'^ 

e(9 - iifi) ~ V2E(9)-- [ d9' ( -(^-^ + -^r^)e^ - ^^^e"^ ) L(9') 
^ ^' n J \^2^sinh^' cosh6''^ sinh26'' / ^ ^ 



TT J 2 I cosh 26*' cosh 26*' ' ^ ^ 



(F.22) 



Notice that at the order which matters, the expansions of e and e are simply propor- 
tional: e_i ~ -\/2e_i and ei ~ V^ei. The proportionality implies that Aperiods + ^/ree = 
|e_iei-\/2(e*'^/^ — e~*'^/^) = — e_iei We can read off the relevant pieces of the expansion 

e_i = E_i - / d9'———L(9') - — / d9'—-—L(9') (F.23) 
V^ttJ sinh2^' ^ ^ 27r7 cosh26'' ^ ' ^ ' 

ei = Ei + / d9' , L(9') + — / d9'— — -L(9') (F.24) 
y/2n J sinh2^' ^ ' J cosh2^' ^ ^ ^ ^ 



where 



E-r = --^ log(6i) + '-e-'-/Hog{l/u2 - 1) (F.25) 
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and 



Ei = ^ log(6i) - -e-/^ log(l/«2 - 1) (F.26) 

Actually, the most useful expression for ^periods + ^ free is the average of the expression 
derived from the large 9 asymptotics and the expression derived from the small 6 asymp- 
totics. Setting e ~ e^e"^ + e^^e^ + • • • and a similar expression for e, which again equals 
-\/2e at this expansion order, 

e' 1 ~ E_i - / de'4-——L{9') + / d9'———L{e') (F.27) 
~ V^n J sinh2^' ^ ^ 27r 7 cosh2^' ^ ^ ^ ^ 

e[ r^Ei + \ de'^——L{e') - — / de'———L{e') (F.28) 

^/2tt J sinh26'' ^ ^ 2tt j cosh26'' ^ ^ ^ ^ 

Putting all together, we can organize the final result in terms with two, one or zero 

contour integrals: eie_i = A2 + Ai+ Aq. 

Ao = {-i/V2\og{bi) +i/2e-'^/^log{l/u2 - l))(z/^log(6i) - i/2e'^/Hog{l/u2 - 1)) 

(F.29) 

simplifies to 

Ao^\ log' h + ^ log'ih - 1/bi) (F.30) 

In the next term, we can average the two available expressions and split again Ai = 
^1,1+^1,2 

1 /" / \ cosh 20' 

Ai,i = -^jde' (-ix/21og(6i)cosh^' +ilog(l/«2 - l))cosh(^' -i7r/4)j ^-^L{e') 

(F.31) 

1 /" / \ sinh 2G' 

^1,2 = ^ j dO' [-^^\og{h^) cosW + i log(l/«2 - 1)) cosh(^' - m/4) j — ^L(^') 

(F.32) 

Next we would like to replace the term in parenthesis with the 9 derivatives of e 
and e respectively, to allow for a useful integration by parts. The difference between the 
terms in parenthesis and the 9 derivatives will generate new terms to be added to A2: 
M,i = ^'1,1 + A2,i + A2,2 and ^1,3 = 3 + ^2,3 + ^2,4- 

1 r msh 29' ~ 
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' V2tv J V2tv J \cosh{e ~ 9') smh29' J ^ sinh2^ ^ ^ 

1 f f jnf^ f sinh26'smh(6'-6'') \ cosh 26* ~ ,^ 

^'^ 27r 7 TT y ^ V cosh 2(^-^0 sinh 2^7 ^ cosh20 ^ ^ 

^'^ 27r y 27r y ^ Vcosh(6' - 6'') cosh2^7 ^ ^ cosh 2^ ^ ' 

On the other hand, we can decompose A2 = ^2,5 + ^2,6 + ^2,7, averaging the two 
expressions and symmetrizing 

v^Try y27ry sinh26'sinh26'' V ; v; V ; 

. _^ /-,,,^inh2(» + ^-)cosh(^-»') 

TT y 27r y cosh2^'sinh2^ v ; v ; v ; 

27r y 27r y cosh 2^ cosh 2^' ^ ^ v ; v ; 

Combining pieces together vast simphfications occur 

^ ^^'^ = 7^1'"^ I '' ^inh2».inh2^-co.h(»-.0 "'^''""" 

+ ^« + ^« - f / ^4 / '' ^o.h 2.- jril;!^.- 2.0 

On the other hand, the pieces to be integrated by parts give 

Here Lz2 indicates Mathematica PolyLog[2,x]. 

In simimary, the final result for Aperiods + A free is given by ( [OOD plus the critical 
value of the Yang- Yang functional which is the sum of( [l:'\42| )-( [b\46[ ). This is the Yang- Yang 
functional for the integral equations in (|F.17| ) ( [F.18|) , written in terms of the new variable 



du = — ^-^^^^20 "^hich makes the kernels symmetric. A shift in the integration contour in 
( ^^461) takes it to the form in ( F^Sl) . 
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F.2. Easy pieces of the hexagon 

Here we focus on the terms of the strong couphng answer for the hexagon which 
contain no integrals. In particular we focus on such terms for the reminder function. 
These terms come from (|1^\3C1| ) plus ABDS-Uke — Abds-, given in [^. This giveJii 

Reasy = -\ (log' + \og\h^ - l/h^)) +\Y. (^^s' + 2LZ2(1 - U,)) (F.47) 

i 

we have ui = -f^, U2 = us = We can write everything in terms of 61, 63 and 

fj, = — 2cos^ by using 62 = ^^'^II^i^Jl'^i^ ■ The limit we are interested in is then bi and 63 
of the same order and very large with (p fixed. One can explicitly check that the above 
combination is finite in such limit. One question we would like to answer is whether we 
could get problematic terms of the form U2 log'^ U2 in the small U2 expansion. Such terms 
could only come from the following piece in the above answer (remember that ui +us = 1 
in the limit.): 

Reasy-part = (log' h + log'(63 " l/h)) + ^ (log' ^2 + 2LZ2(1 - ^2)) (F.48) 

Note that this can be entirely written in terms of bi and 63. In order to show that this 
does not contain the kind of terms mentioned above we use: Li2{l — U2) — —Li2{u2) — 
logW2 log(l — U2) + Furthermore, we call = then we obtain 

Reasy-part = ^ (log^(l " ^2) - 2 log(l - U2) log Z + 2 log^ Z - Li2U2) (F.49) 

From this is clear that in the small U2 expansion we will not get problematic terms, just 
simple power series. Furthermore, in this form it is very easy to expand. 

F.3. Expanding the hexagon at strong coupling 

We choose the following six points, P = {Z^i, Zq, Zi, Z2, Z3, Z4), as 

Pi = (1, e-^\ e-2^ 1, 0, 0) , Ps = (e-2-, 1, -1, -e-^\ 0, 0) , 

P2 = (sinh(r + a), sinh(r — a), — cosh(r — cr), cosh(r + cr), —iV2cos(j), i\/2sm(j)) 

Pe = {l, 1, 1, 1, 0, 0) , P4 = (1, 1, -1, -1, , 0, 0) , P5 = (0, 0, 1, -1, iV2, 0) 

(F.50) 

Recall the relative sign between the remainder function and the area R ~ logVF oc —Area. 
We are ignoring a factor of \/A/(27r). 
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We then define Ui 
These are 



i + l, 1+4^*1 + 2,1 + 5 



62 = 2 (cosh r cosher — COS ^)/ sinh^ r , 6i = e'^coshr, 63 = 6 "^coshr (F.51) 



This imphes that fx defined via fx + = 6162^3 — bi — b2 — bs = —2 cos (f), is n = — e"^. 

The remainder function is = — ^(Area). We wiU suppress the factor of ^ but we 
will take into account the minus sign. Let us now expand the various terms to the desired 
order. Inserting ([F.51|) and (|4.1| ) into (|F.47| ) and expanding we get 



-Reasy ~ 



COS (f)e [cosh a log [2 cosh a] — a sinh a) + 



+ e 
9{^) 



-2f 



log [2 cosher] — a 



+ cos 2(j)g{a) 



(F.52) 



- - cosh 2a log [2 cosh a] + -a sinh 2a -\ — 
4 4 8 



where e""^ = e~^/2. There is a constant term that we neglected. Note that individual 
terms in ( [F.47 ) diverge, but the combination is finite, which is what we expected for the 
remainder function. 

The leading order expansion of the Yang- Yang functional comes only from the di- 
logarithm piece. Expanding the di-logarithm in ( [t'\45| ) ( |1:'\46| ) to first order we get 



Ryy 



YY,r ~ - 



d9 



sinh^ 26 cosh^ 26 



(F.53) 



E{6) = V2[f cosh ^-zcT sinh ^] 



We first used the leading order expressions for e{6 — i(fi) in ( [l''.17| ),( [F'.18D , to get the first 
line. The second line is the approximate expression for (|FTTB]) using the leading order 
expressions for the cross ratios 



- log 2/2 = log[^^ 1] ~ - log U2 ~ -2r 

log 61 = -- logyi^s ~ r + cr 



(F.54) 



Of course, in deriving (|F.52| ) it was important to go beyond this leading order expression. 
( |t'\52| ) and ( |1'\53| ) are the expansions quoted in the main text (|4.2| ), up to the trivial 
relabeling f — )■ r. 
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